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Resource heterogeneity can facilitate cooperation 
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Although social structure is known to promote cooperation, by locally exposing selfish agents 
to their own deeds, studies to date assumed that all agents have access to the same level of 
resources. This is clearly unrealistic. Here we find that cooperation can be maintained when 
some agents have access to more resources than others. Cooperation can then emerge even 
in populations in which the temptation to defect is so strong that players would act fully 
selfishly if their resources were distributed uniformly. Resource heterogeneity can thus be 
crucial for the emergence and maintenance of cooperation. We also show that resource 
heterogeneity can hinder cooperation once the temptation to defect is significantly lowered. 
In all cases, the level of cooperation can be maximized by managing resource heterogeneity. 
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The emergence of cooperation in human and animal 
societies continues to intrigue evolutionary biologists 1 
and sociologists. The fundamental social dilemma 
underlying the tragedy of the commons 2 is described by the 
Public Good Game (PGG) 2 . In the PGG, each of the n players has 
a budget Ri of resources, from which he or she can contribute to a 
common pool that promotes a public good. The pooled resources 
are multiplied by r (l<r<«) and are divided evenly among all 
players (r thus measures the gains from promoting the public 
good). When the group consists only of two players (n — 2) and a 
player can either cooperate by investing all personal resources 
into the common pool, or defect by investing none, this yields the 
following payoffs, denoted in analogy with the Prisoner's 
Dilemma (PD) 3 : T=R 1 + R 1 r/2, R = Rin P = R X and S = R 1 r/2. 
After rescaling (by setting the baseline and unit of payoffs to R\rl2) 
and substituting r — 2/b y this gives T=b, R=l, P—b—l and 
S — 0 y which is a special case (sometimes called 'additive') of the 
payoff structure of the PD 4 . Here the temptation to defect is 
measured by b — I; even though the players are always best-off 
collectively when they cooperate, if b > 1 then an individual 
player will obtain a higher payoff by not contributing, no matter 
what the other players do. Therefore, defection is the 
evolutionarily stable strategy in this game. Although the PGG 
and PD illustrate that cooperators are prone to exploitation and 
that adaptation will favour free riders, ample examples of 
altruistic behaviour exist in human societies 5 and among other 
organisms 6,7 . 

To overcome this incongruity, the simple PGG model can be 
extended and made more realistic in a variety of ways. For 
example, repeated interactions and conditional strategies 8 , 
voluntary participation 9 , costly punishment of defectors 10 , confine- 
ment to groups and elevated mortality of defectors under adverse 
environmental conditions 12 ' 13 , have all been shown to facilitate 
the spread and maintenance of cooperators. Furthermore, em- 
bedding the cooperation game in ecological interactions may 
enable individuals to escape the dilemma 14 ' 15 . Another notable 
extension situates players in physical space or on a social network; 
in such structured populations, cooperators can spread and persist 
if the temptation to defect is not too large 4 ' 16-18 . The reason is that 
local interactions cause defectors to suffer from their own deeds; 
they can locally spread by exploiting cooperators, but then end up 
mostly interacting with other defectors, which drastically reduces 
the advantages of defection. Conversely, cooperators will mostly 
interact with other cooperators, which increase their payoff. 
Recently, interaction topologies other than regular graphs have 
been explored 19 , and it was shown that heterogeneous degree 
distributions facilitate the spread of cooperators 20 . 

Interestingly, previous analyses of the PGG always assumed 
that all players in a locale have access to the same level of 
resources. Therefore, the implications for cooperation of hetero- 
geneity in the resource access, or wealth, of players have so far 
remained unexplored. For example, articles in a recent special 
issue in the Journal of Theoretical Biology devoted to cooperation 
studies do not touch this subject 21 . Although the effect of 
diversity in social structure 22 and of temporally variable resource 
availability for an entire population are discussed, the 
fundamental question of how the diversity of resource 
availabilities within a population affects the level of cooperation 
has not yet been addressed. As uneven distributions of resources 
or wealth are ubiquitous, and as resource heterogeneity has 
already been shown to have an effect in a model of mutualism 24 
between two different species, there is every reason to examine 
how evolutionarily stable levels of cooperation are influenced by 
resource heterogeneity. 

The notion of pooling and dividing resources in the PGG 
naturally generalizes to heterogeneous environments. We assume 
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that each player has access to an amount of either Ri or R 2 of 
resources (with R 2 >Ri). Depending on a player's strategy, these 
resources can be either pooled (cooperator) or retained (defector). 
As before, the pooled resources are multiplied by r, with each 
participating player receiving half of the result. Only two 
parameters suffice to characterize this game (for details, see the 
derivation of the payoff matrix in the Methods section): b — 2lr 
(\<b<2) and c — (R 2 — Ri)/Ri, with b—\ measuring, as before, 
the temptation to defect and c measuring the wealth contrast 
between players on rich and poor sites. We consider a population 
of interacting and competing players that is spatially distributed 
on a square lattice, with rich and poor sites allocated at random 
and the fraction of rich sites given by p\ players interact with their 
four nearest neighbours (other interaction structures are also 
explored). 

We show that cooperation can be facilitated when some agents 
have access to more resources than others. Cooperation can then 
emerge even in populations in which the temptation to defect is 
so strong that agents would act selfishly if their resource access 
were uniform. However, resource heterogeneity can also hinder 
cooperation once the temptation to defect is significantly lowered. 
Our results thus demonstrate how managing resource hetero- 
geneity can be crucial for safeguarding the emergence and 
maintenance of cooperation. 

Results 

At high temptation, heterogeneity facilitates cooperation. Our 

results show that given the same high level of temptation to defect 
(corresponding to relatively low gains from promoting the public 
good), cooperators persist in the heterogeneous environment 
(Fig. lb) even if they would go extinct in a homogeneous 
environment (b —\ >0.25) (ref. 25, Fig. la). In fact, a little bit of 
heterogeneity goes a long way towards stabilizing cooperation 
(Fig. 2a,c,d). 

Although a quantitative analytical theory of this spatial game 
seems beyond reach, some important qualitative insights can 
nevertheless be secured through a combination of numerical and 
analytical examinations. For this purpose, we quantify the most 
important departures between an observed spatial pattern and the 
corresponding mean-field pattern (with the latter defined by the 
random redistribution of site qualities and strategies). Important 
departures are those that are not only large, but also prevalent; we 
measure this by what we call the prevalent- deviation score (PDS) 
f\og 10 (f/f) calculated from a configuration's observed frequency 
/ and mean-field frequency /. By observing the frequencies of 
every configuration of strategies and site qualities involving two 
to five players, and by determining the resultant PDS values 
(Supplementary Note 1), we reveal the common characteristics of 
configurations that are prevalent and strongly deviate from mean- 
field expectations (Fig. 3). 

At high temptation to defect, a rich cooperator is always found 
at the core of a cooperating cluster (Figs lb and 3a,b). In the 
star-shaped configuration that has the highest PDS (Fig. 3b), if 
c/4>b— 1 the poor-site cooperators have higher payoffs than 
any surrounding poor- site defector, whereas in a homogeneous 
environment (c = 0) this condition is never met for b — 1 > 0. 
Hence, cooperators on rich sites enhance the payoffs of 
neighbouring cooperators on poor sites, and thus allow coopera- 
tion to spread in a heterogeneous environment. Our PDS analysis 
in Fig. 3a,b thus suggests that cooperators holding on to small 
rich islands are the main mechanism through which resource 
heterogeneity facilitates the maintenance of cooperation. This 
also explains why the highest facilitation of cooperative behaviour 
occurs at low values of p (Fig. 2a,d): under those conditions, 
rich islands are small, and, thus, are likely to harbour only 
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Figure 1 | Snapshots of the spatial configurations of cooperators and 
defectors. A 30 x 30 snapshot of spatial configurations of cooperators 
(blue) and defectors (red) of the 100 x 100 field in the spatial two-person 
PGG with (b,d) and without (a,c) wealth inequality, with rich sites indicated 
with brighter colours. Site qualities and initial strategies are set 
independently and at random, with (a,c) a p = 0 and (b,d) p = 0.2 fraction 
of rich sites and a 0.5 fraction of cooperators. Portraying the homogeneous 
environment as uniformly rich (p = 1) or poor (p = 0) is equivalent and can 
be used interchangeably. The initial distribution of strategies quickly 
changes within 10 6 steps (middle column), but equilibrium is reached only 
after about 5 x 10 7 steps (b-d, right column). When the temptation to 
defect is high (b — 1 = 0.7; a,b), cooperation cannot be maintained in the 
homogeneous environment (a), whereas it persist in the heterogeneous 
environment (b). Note that clusters of cooperators always harbour a rich 
cooperator at their core (b). When the temptation to defect is low 
(b — 1 = 0.1, c,d), defectors cannot survive in the homogeneous 
environment (c), but they cannot be ousted from the population in the 
heterogeneous environment (d). Note that rich defectors at the core of 
clusters make it possible for defection to be maintained (d). Other 
parameters: c = 10 and m = 10. 

cooperators, whereas the large, rich islands occurring when p is 
high are likely to harbour both cooperators and defectors, and 
consequently are easily taken over by defectors. 



At low temptation, heterogeneity hinders cooperation. When 
the temptation to defect is low (corresponding to relatively high 
gains from promoting the public good), wealth inequality can 
lower cooperation levels (Fig. 2e) compared with a homogeneous 



environment. A low, but not vanishing fraction p of rich sites 
(Fig. 2b,e) and a high contrast c exacerbate this reduction in 
cooperation (Fig. 2e,f). Again, the five-site configuration with the 
highest PDS is star-shaped but now comprise poor-site defectors, 
with a rich-site defector at the centre (Fig. 3d). The rich site 
stabilizes the others if c> \.75/b — 1. Even for low temptations, 
0<b— 1<0.25, this inequality can be satisfied already at low 
contrast, for c>0A to O0.75. As for high temptation, the 
occupancy of small, rich islands is important (Figs lb and 3c,d); 
now, defectors can cling to these, but are ousted from larger 
ones. Analytical results obtained by investigating the dynamics of 
small clusters 25 corroborate the numerical result, by demon- 
strating how such clusters (Supplementary Figs SI and S2) behave 
differently in homogeneous versus heterogeneous environments 
(Supplementary Note 2). 



Results are robust to alternate model assumptions. The above 
results prove to be very robust. Changing the update rule (con- 
sidering synchronous and/or deterministic updating instead of 
the asynchronous stochastic updating employed for the main 
figures) does not affect our qualitative results (Supplementary 
Fig. S3). With deterministic updating, cooperators, at high 
temptation to defect, achieve slightly higher abundance than with 
stochastic updating (Supplementary Fig. S3a). This is in agree- 
ment with other theoretical investigations of cooperation in 
structured populations 17 . At low temptation to defect, the 
prevalence of defectors is much higher than with stochastic 
updating (Supplementary Fig. S3b). This is because a single 
defector amidst cooperators always has a higher payoff than its 
neighbours, so the defector can always spread. With stochastic 
updating, a single defector can go extinct just by chance. This 
explains why full cooperation cannot be achieved with determi- 
nistic updating (Supplementary Fig. S3b,f), but readily results 
for stochastic updating (Fig. 2e and Supplementary Fig. S3d). 
Introduction of non-overlapping generations (synchronous 
updating) causes no significant quantitative deviation from the 
results obtained for overlapping generations (asynchronous 
updating; compare Fig. 2d,e with Supplementary Fig. S3c,d, and 
Supplementary Fig. S3e,f with Supplementary Fig. S3a,b). 

Considering the Moore neighbourhood (eight neighbours), a 
random regular network or a scale-free network result in the same 
qualitative outcomes (Supplementary Fig. S4). The impact of the 
enlarged (Moore) neighbourhood is small. The level of coopera- 
tion increases in the heterogeneous high-temptation environment, 
when the interactions take place on a random regular graph. In 
other words, the randomization of neighbours further increases 
the effect of heterogeneity, so that for high temptation even higher 
levels of cooperation are achieved (Supplementary Fig. S4c,d). 
In contrast, the heterogeneous degree distribution of the scale-free 
interaction graph favours defection in heterogeneous environ- 
ments (Supplementary Fig. S4e,f). Therefore, the levels of 
cooperation achieved in the heterogeneous high-temptation 
environment are more modest compared with those observed 
on regular lattices, and the drop in the levels of cooperation in the 
heterogeneous low- temptation environment is larger. 

Moreover, the possibility of erroneous execution of strategies 
does not affect the qualitative outcomes (Supplementary Fig. S5). 
We can thus conclude that the noise introduced by such errors 
does not have any qualitative effects. When both strategies can 
err, the quantitative results hardly change for high temptation 
(Supplementary Fig. S5a), whereas the negative effects of 
environmental heterogeneity are enhanced when the temptation 
to defect is low (Supplementary Fig. S5b). At intermediate 
frequencies of rich sites, erring defectors can reap a high payoff by 
mistakenly cooperating if there are enough cooperators around. 
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Figure 2 | Fraction of cooperators in heterogeneous environments. Equilibrium fractions of cooperators are shown as a function of the fraction p 
of rich sites and the contrast c between rich and poor sites. Panels a and b show mean levels of cooperation as a function of p of rich sites at c = 9. 
Panels c and f show mean levels of cooperation as a function of the contrast c between rich and poor sites at p = 0.1. Panels a-c and f show transect of the 
parameter space explored in panels d and e, where the location of these transects is indicated by grey lines. Panels d and e depict the fraction of 
cooperators through colour scales from 0 (red) to 1 (blue), based on 10 (a,c,d) or 40 (b,e,f) replicate model runs at each point, with increments of 
0.1 in p and of 1 in c. White lines represent iso-wealth curves; along these, the average resource level of sites, pft 2 + 0 -p)#i, remains constant. Initial 
strategies are set at random, with a 0.5 fraction of cooperators. Results are recorded after 5 x 10 8 steps. Other parameters: m = 10 and b-1 = 0.7 
(high temptation: a,c,d) or b-1 = 0.1 (low temptation; b,e,f). 
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Figure 3 | Configurations with the highest prevalent-deviation score in heterogenous environment. First ten highest-ranking configurations are shown 
for (a,b) two, (c,d) three, (e,f) four and (g,h) five players. The model is first run for 5 x 10 7 steps, after which configurations are recorded every 5 x 10 6 
steps for a total of 20 recordings. Frequencies belonging to mirror-symmetric and rotationally symmetric configurations are averaged and reported for one 
configuration arbitrarily chosen from the symmetry class. Colour coding as in Fig. 1. Other parameters: m = 10 and b — 1 = 0.7, p = 0.1, c = 20 (high 
temptation; a,c,e,g) or b -1 = 0.1, p = 0.3, c = 6.67 (low temptation, b,d,f,h). 



Such players, however, would still spread their defecting strategy. 
When only cooperators can err, the abundance of defectors is 
roughly the same for low temptation as in the absence of 
errors (Supplementary Fig. S5d). Cooperators, however, are 
generally better off if they can occasionally reap the benefit of 



defecting, while still propagating their cooperative behaviour 
(Supplementary Fig. S5c). 

Considering continuous resource levels does not even dis- 
cernibly change our quantitative results (Supplementary Fig. S6). 
The noise introduced by continuous resource levels has negligible 
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effect on the outcome, showing that our assumption of only 
two resource levels does not affect the generality of the reported 
results. 

As a further test of robustness, we analyse the heterogeneous 
PGG for n — 5 players, again finding qualitatively equivalent 
behaviour (Supplementary Fig. S7). Moreover, when we replace 
the PD with the Snowdrift Game 26 ' 27 , an important alternative 
for modelling cooperation, we find that resource heterogeneity 
always favours cooperation (Supplementary Fig. S8). 

In summary, we have shown that our results are robust 
to a wide variety of changes in key features of our model 
(Supplementary Table SI). 

Cooperation can emerge from scratch. Our results show that 
cooperation can be maintained in a heterogeneous environment 
when this is impossible in the corresponding homogeneous 
environment. So far, however, we have assumed that players 
either always cooperate or always defect. In an attempt to 
encompass more real-world complexity, we now extend our 
model by allowing gradual evolution in the individual-level 
probability of players to cooperate. We can thus show that 
cooperation is not only maintained, but can also emerge from 
scratch through gradual evolution in heterogeneous environ- 
ments initially comprising only defectors (Fig. 4a). In this mixed- 
strategy model of gradual evolution, the population-level averages 
of the probability to cooperate turn out to be closely similar to the 
frequencies of cooperators observed in the simpler, pure- strategy 
model (compare Fig. 2d,e with Fig. 4a,e, respectively). 

These results can be understood by appreciating that the 
cooperative act is an investment towards the objective of being 



surrounded by cooperators. Players on rich sites can maintain 
this cooperative neighbourhood by being cooperative themselves, 
up to the point at which the resultant advantage is balanced by 
the danger of facilitating invasion by a less cooperative player 
residing on an adjacent rich site. 

In a further extension, we analyse the evolution of conditional 
strategies, with each player now being characterized by two 
continuous variables: the probability to cooperate if the player is 
on a rich site, and the probability to cooperate if the player is on a 
poor site. Allowing conditional strategies to evolve results in 
amplified effects of heterogeneity. With conditional strategies, 
only a little contrast between rich and poor sites (c= 1) already 
suffices for the population to reap the full benefit of wealth 
inequality through increased cooperation (Fig. 4b). Moreover, the 
range of fraction p of rich sites over which heterogeneity enables 
cooperation is considerably broadened relative to the uncondi- 
tional case (Fig. 4b). Remarkably, players evolve a higher 
probability to cooperate when situated on poor sites than when 
residing on rich sites (compare Fig. 3c,d), with the probability to 
cooperate on rich sites increasing, as rich sites become scarce. 
Note that when there are very few rich sites, the players on these 
sites are more inclined to cooperate and the resultant average level 
of cooperation is mainly due to their contributions (compare Fig. 
4c,d). At low levels of temptation, populations cease to exhibit full 
cooperation when poor sites abound, and this departure grows as 
wealth inequality increases (Fig. 4f). In this case, players on rich 
sites are reluctant to cooperate when their numbers are low, as 
exploitation of their altruistic behaviour leads to elevated 
defection. Consequently, the probability of players on rich sites 
to cooperate steadily increases with their abundance (Fig. 4g). 
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Figure 4 | Evolved propensities to cooperate in heterogeneous environment. Average equilibrium probabilities to cooperate are shown as a function of 
the fraction p of rich sites and the contrast c between rich and poor sites. Panels a and e show the evolutionary outcomes for unconditional strategies, 
whereas panels b-d and f-h show those outcomes for conditional strategies. In the latter case, the probabilities to cooperate when situated on a rich 
or poor site evolve independently. Average equilibrium probabilities to cooperate are shown in panels b and f. The mean equilibrium probabilities of 
cooperation of players on rich sites are shown on panels c and g, and the same for players on poor sites in panels d and h. All panels depict the fraction of 
cooperators through colour scales from 0 (red) to 1 (blue), based on five replicate model runs at each point, with increments of 0.1 in p and of 1 in c. Initial 
strategies were set to 0, implying defection on both rich and poor sites. As in Fig. 2, white lines represent iso-wealth curves. Results are recorded after 
2 x 10 9 steps. Other parameters: m = 10 and b — 1 = 0.7 (high temptation; a,c,d) or b — 1 = 0.1 (low temptation; b,e,f). 
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In general, players on poor sites tend to be more cooperative than 
players on rich sites. This is because a group of poor-site players can 
increase their payoffs more if they become more cooperative, 
whereas a group of rich-site players can maintain their competitive 
superiority even with a more modest cooper ativeness. 

Discussion 

Resource heterogeneity allows cooperation to persist even when 
the temptation to defect is so high that full defection is observed 
in the homogeneous case. This key finding of our study is broadly 
related to an earlier research, demonstrating that asymmetries 
among players can facilitate the spread of cooperation. Maynard 
Smith , for instance, showed that role asymmetry through the 
introduction of the 'bourgeois' strategy can lessen aggression in 
the Hawk-Dove Game. Moreover, variation in the total number 
of interaction links individuals possess on heterogeneous 
interaction graphs results in elevated levels of cooperation 20 ' 29 . 
This is because some individuals can then reap benefits from 
many interactions, resulting in inequalities in the maximal 
benefits of players 30 . Similarly, the existence of role models 
(individuals whose successes are copied more often than those of 
others) can stabilize cooperation 31 ' , but see Kim et al. 33 Thus, 
differences among individuals appear to be an important, and as 
yet insufficiently studied, determinant of the cooperation level a 
society can attain. 

Interpreting our results in a wider context, it is tempting to ask 
whether they might contribute an explanatory facet to the 
scientific understanding of more or less egalitarian social 
norms 34 . Our results suggest that the level of cooperation can 
be managed by changing the distribution of wealth (moving along 
the iso -wealth curves in Figs 2 and 3). Although differences 
among individuals are pervasive in the living world and in human 
societies, social norms can either counter or exacerbate them. 
Specifically, in situations in which the temptation to defect is high 
(as gains from promoting the public good are low), it will always 
be possible to increase the levels of cooperation by increasing 
wealth inequality (Figs 2d and 3a,b). On the other hand, 
in situations in which gains from promoting the public good 
are high (and the temptation to defect thus is low), wealth 
inequalities lead to decreased levels of cooperation (Figs 2e and 
3e,f). In the latter case, the evolution of inequality aversion, 
promoting an egalitarian society, can considerably increase levels 
of cooperation by decreasing wealth inequalities, ultimately 
leading to a situation in which all members of a population 
receive the same share of resources. 

Egalitarian motives in humans are widespread and are 
exhibited from a very early age 35-37 . Fehr and Schmidt 38 have 
thus suggested equity as a key factor in strengthening human 
cooperation and in explaining the prevalence of cooperation in 
our species 34 ' 39 . The case for animals having inequality aversion 
is less strong. Some studies have documented strong reactions by 
animals towards what they perceive as inequity 40-43 ; yet, the 
larger shares accrued by socially dominant individuals are seldom 
questioned by other individuals 41 and the resultant social 
disapproval is even less strong towards members of the same 
social group 42 . On the other hand, some studies found that great 
apes are true rational maximizers when it comes to their feeding 
options 44,45 and no inequality aversion could be found among 
them. If we account for the facts that great apes and humans both 
exhibit cooperative behaviour 46 , but human cooperation is much 
more pervasive 47 and efficient at generating public goods, we can 
arrive at the conclusion that (early) humans and great apes are 
positioned on opposite sides of the threshold in the efficiency of 
public-good generation; below this threshold heterogeneity and 
inequity are beneficial, whereas above this threshold homogeneity 



and equity should be the norm (Fig. 2d,e). Examples for such 
disparate efficiencies abound, for example, chimpanzees, can 
cooperatively hunt only small monkeys 48 , whereas human 
hunters can take down large preys 49 ; likewise, chimpanzees 
physically attack only lonely individuals from other bands, 
whereas humans regularly engage in warfare 50 . This suggests that 
cooperation was a necessity in the evolutionary past of our 
species, and that the efficiency of our cooperation has co- evolved 
with our egalitarian motives. 

Although humans appear to have been fully egalitarian at the 
hunter-gatherer stage of cultural evolution 51 , and although 
contemporary hunter-gatherer societies still exhibit low levels of 
wealth inequality 52 , industrial and post-industrial societies have 
varied levels of wealth inequality 53 and also show variation in the 
level of observable cooperative tendencies 54-56 . This suggests that 
our results might help explain the re- emergence of inequality 
during the Holocene . As human societies became larger, and 
mechanisms limited to enabling the maintenance of cooperation in 
small-scale societies (like peer reputation and peer punishment 55 ) 
thus lost their effectiveness, institutions emerged that allowed 
enforcing cooperative norms at a cost 58 ' 59 . Maintaining the same 
level of cooperation at the lower efficiency implied by the involved 
cost required elevated wealth inequality, which indeed prevailed for 
hundreds of years 57 . Today we hopefully are in the middle of a new 
transition in the cultural evolution of our species, back to more 
egalitarian societies. Indeed, modern Western societies exhibit high 
levels of cooperativeness 56 coupled with low wealth inequalities 
comparable to those of the hunter-gatherers 52 . Our model-based 
analyses suggest that this transition could have come about by an 
increased effectiveness of public-good generation. Egalitarianism, a 
unique human trait, can once more hold sway in our societies. 

Methods 

Payoff matrix. The payoff matrix (Table 1) is obtained in four steps. First, players 
on a poor site can invest R lt whereas players on a rich site can invest R 2 ; coop- 
erators invest, whereas defectors do not. The pooled investments are multiplied by 
r and then evenly distributed among the two players. Second, payoffs are nor- 
malized by adjusting the payoff baseline so that a cooperator on a poor site playing 
against a defector on a poor site receives a payoff of 0. Third, payoffs are nor- 
malized by adjusting the payoff unit, so that two cooperators on poor sites playing 
against each other will each receive a payoff of 1. These two normalizations have no 
bearing on game dynamics. Fourth, the definitions b = 2/r and c = (R 2 — R\)/Ri are 
applied to arrive at the payoff matrix shown below. Note that switching coop- 
eration to defection in this matrix always increases a player's payoff by b — 1 on a 
poor site and by (c + l)(b — 1) on a rich site, irrespective of the partner's strategy or 
site quality. Note also that the site quality of partners affects a player's payoff only 
when their partners cooperate, but not when they defect. 

Game dynamics. Our spatial game dynamics are analogous to those analysed, for 
example, by Nowak et al. 4 Individual players reside on sites of a 100 x 100 square 
lattice with periodic boundary conditions. Updating is asynchronous, by randomly 
choosing a focal site to be updated in each step. The player on this focal site then 
plays the game described in the main text with all four nearest neighbours, and 



Table 1 | Payoff matrix of the two-person Public Good Game 
with wealth inequality. 


First player 




Second player 


Rich site 


Poor site 


Cooperate 


Defect 


Cooperate Defect 


Rich site 








Cooperate 


2c + 1 


c 


c + 1 c 


Defect 


bc + b + c 


bc + b-1 


be + b bc + b-1 


Poor site 








Cooperate 


c+1 


0 


1 0 


Defect 


b + c 


b-1 


b b-1 



6 



NATURE COMMUNICATIONS | 4:2453 | DOI: 10.1038/ncomms3453 | www.nature.com/naturecommunications 
© 2013 Macmillan Publishers Limited. All rights reserved. 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms3453 



ARTICLE 



these neighbours in turn play it with their four nearest neighbours. The focal site is 
held by the focal player, or is taken over by its neighbours, with probabilities 
proportional to their payoffs raised to the power of m (resulting in what is known 
as proportional or imitation updating 27 ), which thus characterizes the degree of 
stochasticity in the contest for a site 4 . When a site's occupancy changes, the payoffs 
of the focal player and its neighbours are updated. 

Evolutionary dynamics. In our evolutionary analysis, the game dynamics are as 
described above, except that the lattice size is reduced to 64 x 64 sites to keep 
computation time manageable. Each player's strategy is given by his or her 
probability to cooperate in an interaction. If a site's occupancy changes, there is a 
5% chance that the focal player's strategy is mutated. The mutated strategy is 
chosen randomly from the interval [x — 0.05, x + 0.05], where x is the strategy 
before mutation, and any mutated strategy outside the interval [0,1] is reset to the 
corresponding boundary. Evolutionary dynamics start with every player being a 
defector, implying a vanishing probability to cooperate. We also allow strategies to 
be conditional on the quality of the occupied site, resulting in two cooperation 
probabilities that can evolve independently. Upon mutation, one of these two 
probabilities is chosen at random and is changed as described above. 

Robustness checks. Many alternate model assumptions are examined. With 
asynchronous deterministic updating, the player with the highest payoff always 
wins the local competition, which is equivalent to m — »■ oo . This update rule is 
called 'best takes over' 26 . With synchronous stochastic updating, all players' payoffs 
and strategies are updated at the same time. This corresponds to the dynamics of a 
population with non-overlapping generations. In practice, synchronous updating is 
carried out by determining the outcome of competition based on the state of a site 
at time t, and continuing with the winners of all corresponding local competitions 
at time t+1. With synchronous deterministic updating, the player with the highest 
payoff always wins the local competition and all players' payoffs and strategies are 
updated at the same time. 

In the model variant in which the interaction structure is described by the Moore 
neighbourhood, each player plays with its eight neighbours, which include the four 
orthogonal and the four diagonal neighbours. In the case when interactions take 
place on a random regular graph, the neighbourhood relations are fully randomized 
at the beginning of each model run. This means that each player has k = 4 
neighbours, but these players can be located anywhere on the original lattice. The 
scale-free graphs are generated by the Albert-Barabasi method 60 at the beginning of 
each model run. In this case, on average, each player has k = 4 neighbours, but the 
actual numbers of neighbours follow a power-law distribution. Payoffs are scaled to 
the actual number of neighbours before comparing them in competition. 

Erroneous strategy execution is modelled as follows: in each individual game, 
each player has a small probability (e = 0.01) accidentally to play the opposite 
strategy (that is, cooperators defect, whereas defectors cooperate). Of course, such 
an event only affects strategy execution at one moment but never lastingly changes 
the strategy itself. As a further test, we examine what happens if defectors would 
not, by mistake, behave as cooperators. We, therefore, also consider a setup in 
which cooperators can erroneously defect, whereas defectors never cooperate. As 
players in our model do not base their behaviours on the behaviours perceived in 
their neighbours, erroneous perception need not be considered. 

Continuous levels of resources are incorporated by considering a continuous 
bimodal distribution of resources instead of just two discrete levels of resources. 
Specifically, resources are either low, Ri + p, or high, R 2 + p, where p is a random 
number drawn from a Gaussian distribution with mean 0 and s.d. o. Payoffs are 
calculated by pooling the available resources from cooperators, dividing by b = 21 r 
(see the payoff matrix section of the Methods), and giving the result to both 
participant of the game. Defectors withhold their resources. 

In the heterogeneous Snowdrift Game, each player can gain a benefit from 
cooperation, b x on rich sites or b 2 on poor sites (with b x >b 2 ), and the cost of 
cooperation is c 0 for any player. In the Snowdrift Game, the cost is halved when 
both players cooperate, and the benefit is reaped when at least one of the players 
cooperates. We use synchronous stochastic updating, as described above. 

The five-player PGG is a generalization of the two-player PGG, in which a 
player and N — 1 of its neighbours constitute a group playing a round of the PGG. 
As in the basic game, players on poor sites can invest R lt whereas players on rich 
sites can invest JR 2 ; cooperators invest, whereas defectors do not. The pooled 
investments are multiplied by r and are then evenly distributed among the N 
players. After applying the definition of contrast, c = (R 2 — Ri)/R\, and normalizing 
payoffs by R x , we obtain the payoff of a cooperator as IT = r(A: poor + (c + l)A; ric h)/N, 
where k poor and k rich are the numbers of cooperators within the group situated on 
poor and rich sites, respectively. Analogously, the payoff of a defector is IT + 1 on a 
poor site and IT + (c + 1) on a rich site. 

References 

1. Maynard Smith, J. & Szathmary, E. The Major Transition in Evolution (W.H. 
Freeman, 1995). 

2. Hardin, G. The tragedy of the commons. Science 162, 1243-1248 (1968). 

3. Axelrod, R. & Hamilton, W. D. The evolution of cooperation. Science 211, 
1390-1396 (1981). 



4. Nowak, M. A., Bonhoeffer, S. & May, R. M. Spatial games and the maintenance 
of cooperation. Proc. Natl Acad. Sci. 91, 4877-4881 (1994). 

5. Camerer, C. F. Behavioural Game Theory: Experiments in Strategic Interaction 
(Princeton University Press, 2003). 

6. Nunn, C. L. & Lewis, R. J. Cooperation and collective action in animal 
behaviour, in Economics in Nature (eds Noe, R., van Hooff, J.A.R.A.M. 
& Hammerstein, P.) 42-66 (Cambridge University Press, 2001). 

7. Dugatkin, L. A. Animal cooperation among unrelated individuals. 
Naturwissenschaften 89, 533-541 (2002). 

8. Nakamaru, M., Matsuda, H. & Iwasa, Y. The evolution of cooperation in a 
lattice- structured population. /. Theor. Biol. 184, 65-81 (1997). 

9. Szabo, G. & Hauert, C. Evolutionary prisoner's dilemma games with voluntary 
participation. Rhys. Rev. E 66, 062903 (2002). 

10. Nakamaru, M. & Iwasa, Y. The evolution of altruism by costly punishment in 
lattice- structured populations: score-dependent viability versus score- 
dependent fertility. Evol. Ecol. Res. 7, 853-870 (2005). 

11. Killingback, T., Bieri, J. & Flatt, T. Evolution in group-structured populations 
can resolve the tragedy of the commons. Proc. R. Soc. Lond. B 273, 1477-1481 

(2006) . 

12. Harms, W. Cooperative boundary populations: the evolution of cooperation on 
mortality risk gradients. /. Theor. Biol. 213, 299-313 (2001). 

13. Kun, A. & Oborny, B. Survival and competition of clonal plant populations in 
spatially and temporally heterogeneous habitats. Com. Ecol. 4, 1-20 (2003). 

14. Mesterton- Gibbons, M. An escape from the 'prisoner's dilemma'. /. Math. Biol. 
29, 251-269 (1991). 

15. Worden, L. & Levin, S. A. Evolutionary escape from the prisoner's dilemma. 
/. Theor. Biol. 245, 411-422 (2007). 

16. Doebeli, M. & Hauert, C. Models of cooperation based on the Prisoner's 
Dilemma and the Snowdrift game. Ecol. Lett. 8, 748-766 (2005). 

17. Nowak, M. A. & Sigmund, K. Games on grids, in The Geometry of Ecological 
Interactions (eds Dieckmann, U., Law, R. & Metz, J. A. J.) 135-150 (Cambridge 
University Press, 2000). 

18. Cavaliere, M., Sedwards, S., Tarnita, C. E., Nowak, M. A. & Csikasz-Nagy, A. 
Prosperity is associated with instability in dynamical networks. /. Theor. Biol. 
299, 126-138 (2012). 

19. Szabo, G. & Fath, G. Evolutionary games on graphs. Phys. Rep. 446, 97-216 

(2007) . 

20. Santos, F. C, Pacheco, J. M. & Lenaerts, T. Evolutionary dynamics of social 
dilemmas in structured heterogeneous populations. Proc. Natl Acad. Sci. 103, 
3490-3494 (2006). 

21. Nowak, M. A. Evolving cooperation. /. Theor. Biol. 299, 1-8 (2012). 

22. Santos, F. C, Pinheiro, F. L., Lenaerts, T. & Pacheco, J. M. The role of diversity 
in the evolution of cooperation. /. Theor. Biol. 299, 88-96 (2012). 

23. Tavoni, A., Schliiter, M. & Levin, S. The survival of the conformist: social 
pressure and renewable resource management. /. Theor. Biol. 299, 152-161 
(2012). 

24. Boza, G. & Scheuring, I. Environmental heterogeneity and the evolution of 
mutualism. Ecol. Compl. 1, 329-339 (2004). 

25. Hauert, C. Fundamental clusters in spatial 2x2 games. Proc. R. Soc. Lond. B 
268, 761-769 (2001). 

26. Hauert, C. & Doebeli, M. Spatial structure often inhibits the evolution of 
cooperation in the snowdrift game. Nature 428, 643-646 (2004). 

27. Kun, A., Boza, G. & Scheuring, I. Asynchronous snowdrift game with 
synergistic effect as a model of cooperation. Behav. Ecol. 17, 633-641 (2006). 

28. Maynard Smith, J. Evolution and the Theory of Games (Cambridge University 
Press, 1982). 

29. Santos, F. C. & Pacheco, J. M. Scale-free networks provide a unifying 
framework for the emergence of cooperation. Phys. Rev. Lett. 95, 098104 
(2005). 

30. Szolnoki, A., Perc, M. & Danku, Z. Towards effective payoffs in the prisoner's 
dilemma game on scale-free networks. Phys. A 387, 2075-2082 (2008). 

31. Wu, Z.-X., Xu, X.-J., Huang, Z.-G., Wang, S.-J. & Wang, Y.-H. Evolutionary 
prisoner's dilemma game with dynamic preferential selection. Phys. Rev. E 74, 
021107 (2006). 

32. Szolnoki, A. & Szabo, G. Cooperation enhanced by inhomogeneous activity of 
teaching for evolutionary Prisoner's Dilemma games. Lett. J. Expl. Frontiers 
Phys. 77, 30004 (2007). 

33. Kim, B. J. et al. Dynamic instabilities induced by asymmetric influence: 
prisoners' dilemma game in small- world networks. Phys. Rev. E 66, 021907 
(2002). 

34. Dawes, C. T., Fowler, J. H., Johnson, T., McElreath, R. & Smirnov, O. 
Egalitarian motives in humans. Nature 446, 794-796 (2007). 

35. Schmidt, M. F. H. & Sommerville, J. A. Fairness expectations and altruistic 
sharing in 15-month-old human infants. PLoS One 6 (2011). 

36. Fehr, E., Bernhard, H. & Rockenbach, B. Egalitarianism in young children. 
Nature 454, 1079-1083 (2008). 

37. Blake, P. R. & McAuliffe, K. T had so much it didn't seem fair': eight-year-olds 
reject two forms of inequity. Cognition 120, 215-224 (2011). 

7 



NATURE COMMUNICATIONS | 4:2453 | DOI: 10.1038/ncomms3453 | www.nature.com/naturecommunications 

© 2013 Macmillan Publishers Limited. All rights reserved. 



ARTICLE 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms3453 



38. Fehr, E. & Schmidt, K. M. A theory of fairness, competition, and cooperation. 
Q. /. Econ. 114, 817-868 (1999). 

39. Fehr, E. & Fischbacher, U. The nature of human altruism. Nature 425, 785-791 
(2003). 

40. Brosnan, S. F. & De Waal, F. B. M. Monkeys reject unequal pay. Nature 425, 
297-299 (2003). 

41. Brosnan, S. F. Nonhuman species' reactions to inequity and their implications 
for fairness. Soc. Justice Res. 19, 153-185 (2006). 

42. Brosnan, S. F., Schiff, H. C. & De Waal, F. B. M. Tolerance for inequity may 
increase with social closeness in chimpanzees. Proc. R. Soc. Lond. B 272, 
253-258 (2005). 

43. Range, F., Horn, L., Viranyi, Z. & Huber, L. The absence of reward induces 
inequity aversion in dogs. Proc. Natl Acad. Sci. 106, 340-345 (2008). 

44. Jensen, K., Hare, B., Call, J. & Tomasello, M. What's in it for me? Self-regard 
precludes altruism and spite in chimpanzees. Proc. R. Soc. Lond. B 273, 
1013-1021 (2006). 

45. Jensen, K., Call, J. & Tomasello, M. Chimpanzees are rational maximizers in an 
ultimatum game. Science 318, 107-109 (2007). 

46. Van Schaik, C. P. & Kappeler, P. M. in Cooperation in Ptrimates and Humans 
(eds Van Schaik, C. P. & Kappeler, P. M.) 3-24 ( Springer- Verlag, 2006). 

47. Jaeggi, A. V., Burkart, J. M. & Van Schaik, C. P. On the psychology of 
cooperation in humans and other primates: combining the natural history and 
experimental evidence of prosociality. Philos. T. R. Soc. B 365, 2723-2735 
(2010). 

48. Boesch, C. Cooperative hunting in chimpanzees. Anim. Behav. 48, 653-667 
(1994). 

49. Diamond, J. The Third Chimpanzee: The Evolution and Future of the Human 
Animal (Harper Perennial, 2006). 

50. Crofoot, M. C. & Wrangham, R. W. in Mind the Gap (eds Kappeler, P. M. & 
Silk, J. B.) 171-195 (Springer Verlag, 2010). 

51. Boehm, C. Hierarchy in the Forest (Harvard University Press, 1999). 

52. Smith, E. A. et al. Wealth transmission and inequality among hunter-gatherers. 
Curr. Anthropol. 51, 19-34 (2010). 

53. Sala-i-Martin, X. The world distribution of income: falling poverty and ... 
convergence, period. Q. /. Econ. 121, 351-397 (2006). 

54. Lamba, S. & Mace, R. Demography and ecology drive variation in cooperation 
across human populations. Proc. Natl Acad. Sci. 108, 14426-14430 (2011). 

55. Gachter, S. & Herrmann, B. Reciprocity, culture and human cooperation: 
previous insights and a new cross-cultural experiment. Philos. T. R. Soc. B 364, 
791-806 (2009). 

56. Gachter, S., Herrmann, B. & Thoni, C. Culture and cooperation. Philos. T. R. 
Soc. B 365, 2651-2661 (2010). 



57. Bowles, S., Smith, E. A. & Mulder, M. B. The emergence and persistence of 
inequality in premodern societies: introduction to the special section. Curr. 
Anthropol. 51, 7-17 (2010). 

58. Baldassarri, D. & Grossman, G. Centralized sanctioning and legitimate 
authority promote cooperation in humans. Proc. Natl Acad. Sci. 108, 
11023-11027 (2011). 

59. Sasaki, T., Brannstrom, A., Dieckmann, U. & Sigmund, K. The take-it-or-leave- 
it option allows small penalties to overcome social dilemmas. Proc. Natl Acad. 
Sci. 109, 1165-1169 (2012). 

60. Albert, R. & Barabasi, A.-L. Statistical mechanics of complex networks. Rev. 
Mod. Phys. 74, 47-97 (2002). 

Acknowledgements 

This work was supported by grants from the Hungarian National Research Fund (OTKA 
D048406, K100299) to A.K. and from the Austrian Science Fund (FWF, M983-N18) to 
A.K. and U.D. A.K. gratefully acknowledges the Bolyai Janos fellowship of the Hungarian 
Academy of Sciences. U.D. gratefully acknowledges additional support by the Sixth 
Framework Program of the European Commission, the European Science Foundation 
(FWF, TECT 1-106 Gil), the Austrian Science Fund, the Austrian Ministry of Science 
and Research, and the Vienna Science and Technology Fund. 

Author contributions 

A.K. instigated the study, which the authors then designed together. A.K. carried out the 
simulations and analysed the resultant data. U.D. developed the PDS method. The 
authors jointly wrote the manuscript and approved its final version. 

Additional information 

Supplementary Information accompanies this paper at http://www.nature.com/ 
naturecommunications 

Competing financial interests: The authors declare no competing financial interests. 

Reprints and permission information is available online at http://npg.nature.com/ 
reprintsandpermissions/ 

How to cite this article: Kun, A. and Dieckmann U. Resource heterogeneity can facilitate 
cooperation. Nat. Commun. 4:2453 doi: 10.1038/ncomms3453 (2013). 



This work is licensed under a Creative Commons Attribution- 
NonCommercial-NoDerivs 3.0 Unported License. To view a copy of 



this license, visit http://creativecommons.Org/licenses/by-nc-nd/3.0/ 



8 



NATURE COMMUNICATIONS | 4:2453 | DOI: 10.1038/ncomms3453 | www.nature.com/naturecommunications 
© 2013 Macmillan Publishers Limited. All rights reserved. 



